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The on-shell self-energy of the homogeneous electron gas in second order of exchange, 
^2x = Re £2x(^F> kp/2), is given by a certain integral. This integral is treated here 
in a similar way as Onsager, Mittag, and Stephen [Ann. Physik (Leipzig) 18, 71 
(1966)] have obtained their famous analytical expression e2 x = 

l ln2 - |^3) (in 

atomic units) for the correlation energy in second order of exchange. Here it is 
shown that the result for the corresponding on-shell self-energy is T,2 X = C2 X - The 
off-shell self-energy Y,2 x (k,uj) correctly yields 2e2 x (the potential component of e2 x ) 
through the Galitskii-Migdal formula. The quantities e2 x and S2 X appear in the 
high-density limit of the Hugenholtz-van Hove (Luttinger-Ward) theorem. 
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I. INTRODUCTION 



Although not present in the Periodic Table the homogeneous electron gas (HEG) is still an 
important model system for electronic structure theory, cf. e.g. [ifl. In its spin-unpolarized 
version the HEG ground state is characterized by only one parameter r s , such that a 
sphere with the radius r s contains on average one electron jgj. It determines the Fermi 
wave number as kp = l/(ar s ) with a = (4/(97r)) 1//3 and it measures simultaneously the 
interaction strength and the density such that high density corresponds to weak interaction 
and hence weak correlation. For recent papers on this limit cf. 0, U, 0]. Naively one 
should expect that in this weak- correlation limit the Coulomb repulsion ar s /r (where 
lengths and energies are measured in units of kp 1 and fcf>, respectively) can be treated 
as perturbation. But in the early theory of the HEG, Heisenberg 6( has shown, that 
ordinary perturbation theory does not apply. With Cq being the energy per particle of the 
ideal Fermi gas and e x being the exchange energy in lowest (1st) order, the total energy 
e = eo + e x + e c defines the correlation energy e c = ei + e?„ + ■ ■ ■ . In 2nd order, there 
is a direct (d) term e2d and an exchange (x) term e2 X , so that e^ = e2d + e2 X . Whereas 
e2x/(oTs) 2 is a pure finite number 6 X , the direct term e2d logarithmically diverges along 
the Fermi surface (i.e. for vanishing transition momenta q): e2d — > lng for q — > 0. This 
failure of perturbation theory has been repaired by Macke j?| with an appropriate partial 
summation of higher-order terms up to infinite order. The result of this (ring-diagram) 
summation for the correlation energy in its weak-correlation limit is e c / (ar s ) 2 = a\nr s + ■ ■ ■ 
with a = (1 — ln2)/7r 2 ~ 0.031091. This has been confirmed later by Gell-Mann and 
Brueckner a] , who in addition to the logarithmic term numerically calculated contributions 
to the next (constant, i.e. not depending on r s ) term b, namely b r and arising from the 
ring-diagram summation and from e2d, respectively: b r w aln^ — 0.001656 ~ —0.057514 
and 6d ~ —0.013586. The total constant term is b = b r + + b x ~ —0.046921, where the 
exchange term b x = | In 2 — ~ +0.0241792 has been analytically calculated by Onsager, 
Mittag, and Stephen in a very tricky way 0. Although the integrand is rather simple, 
the Pauli principle causes complicated boundary conditions, which need a sophisticated 
treatment through several substitutions of the integral variables. 



Here their method is used to calculate the analog term of the self-energy on the energy shell, 



namely £ 2x = Re £ 2x (l, 1/2) [with k measured in units of k-p and u measured in units of fcf]. 
The more general quantity £ 2x (A;,u;) appears (i) in a recent study of the spectral moments 



of the HEG 



Q, (n) 



in the context of self-consistent GW calculations 



iJ3 



. For k = 1, oj 



|l2j ]. and within a 
\i it appears in the 



project 'spectral functions for a hydrogen plasma 

Hugenholtz-van Hove (Luttinger-Ward) identity 14], which relates the chemical potential 
fi to the self-energy Ti(k,cu) according to \x — /in = £(l,/i). The chemical potential follows 
from the total energy according to // = (| — ^f s -£r)e. In the weak-correlation limit r s — > 
the total energy e = e + e x + e c and the chemical potential fi = fi Q + fi x + fi c are given by 



e 



3 

To' 



/io- ^> 



/' x 



e x = 

7T 



3 ar s 

4 7T 

A* c = 



or. 



a In r s + 



[a lnr s + b + 0(r s )], 



--o + 6 J + 0(r s , 



Because of the above mentioned identity the self-energy £(l,/i) must have a corresponding 
behavior. The exchange in lowest order yields S x (l) = — ^* or /i x = S x (l), exactly in 
agreement with the mentioned identity. To obtain also the lo gar ithmic term of £(1, 1/2) the 
ring-diagram summation has to be done for the self-energy [15]. To the next term beyond 
the logarithmic term contributes the 2nd-order exchange self-energy S 2x — Re £2 X (1, 1/2). 
Just this term is calculated here using the tricky method of Onsager, Mittag, and Stephen 
[l| mutatis mutandi. The Feynman diagrams of e 2x and S 2x are given in Figs. 1 and 2. As 



shown in the Appendix, from the Feynman diagram rules it follows S 2x = -&^(X 1 +X 2 ), 
where Ai )2 mean the integrals defined in Eqs. ()A.10j) and (jA.llJ) . They are calculated in the 



following sections. The final results are X\ 



-7T 



ln2-5^ 



X, 



71 



' I2 ln2 _ 2 # 



thus E 2x 



|ln2-|# 

O 4 7T^ 



e 2x , thus /i 2x = ^2 X - This relation appears in the weak- 



correlation limit of the above mentioned Hugenholtz-van Hove theorem 



3,0- 



This paper 



is a contribution to the mathematics of the weakly-correlated (high-density) HEG. 



II. THE INTEGRAL X x 

In Eq. (jA.lOp . new variables q' = (q 1 + g 2 )/2, p' = (q 1 — <? 2 )/2, and s' = e + (q 1 + q 2 )/2 
lead to 



qi = q' + p', (e + Ql f = (s' + p')\ (e + q 2 f = (s> - p')\ 

q 2 = q'- p\ (e + q 1 + q 2 f = (s' + q')\ e 2 = (s' - q'f . 



(2.1) 



Therefore (with Q(x) = Heaviside step function) 

> 3 1 



X 1 = I d 3 q'd 3 p' 



p' 2 {p + q') 2 (p' - q') 2 
x J ^ 5 (( s ' - <i') 2 - i) ! 1 - ( s> + q') 2 M(s' + p'f - i]e[(«' - P ') 2 - 1]. (2.2) 



The next step scales g', p', s' with A = 1 / vl — s /2 according to q = Aq', p = Ap', s = As' 
with the consequences 

1 d 2 s 

i- s ' 2 = — L_, d 2 s'- d 



1 + s 2 ' (1 + s 2 ) 2 ' 

5((s'-q') 2 -l) = (l + s 2 ) 5(q 2 -2sq-l), 

±2sp + p 2 >l, 2sq + q 2 <l, -2sq + q 2 = 1, (2.3) 

from which follow p > 1 and g < 1 (what makes the energy denominator q 2 — p 2 negative) 
and sq < 0, s > a. Thus 

1 oo +1 

/",/", 16 7r f dx 8 q 2 p 2 
X 1 = - dq dp — / — . , ■ ,.- — — - x 



p 2 - q 2 J 2 (p 2 + g 2 ) 2 - (2qpy 

1 -1 

r d 2 s 

x / 5{q 2 -\-2sq) Q{p 2 -l + 2sp) Q(p 2 -l-2sp). (2.4) 

J 1 + s 

s>a, COS(p q <0 

Here and in the following of Sec. II the abbreviations 

1 — q 2 v 2 — 1 q 2 + V 2 

a = —r^ L >0, P = ^-= >0, a = ^-P->l, 

2q 2p 2qp 

- = COS (p p = COS <(s,p), X = COS ip = COS <[ (q, p) , y = cos</? ? = cos<£(q, s) (2.5) 
are used: 

1 oo +1 

f 7 f 7 16 n f dx f sds difia 1 . /a \ ... 

1 —1 s>a, cosip q <0 

t is introduced to replace s after having done the ip q integration. To this purpose the 
relation between t and s is needed. It follows from ip q + ip p = ip or ip p = p — ip q , therefore 
cos ip p = cos((/9 — tp q ) = cos ip cos ip q + sin ip sin ip q , what is in terms oit/s, x, and y (the latter 
quantity equals — a/s because of the delta function): 



- s = x (-^) ± y/T^^l - or (t + ax) 2 = (l-x 2 )(s 2 -a 2 ). 



(2.7) 



This has the consequences 



2 2 
s — a 



(t + ax) 



1 -x 2 



sds 



(t + ax)dt 



(2.8) 

1 + s 2 (t + ax) 2 + (1 + a 2 )(l - x 2 ) ' ^ 
With the help of Eq. ()2.8|) the (p q integration yields a function depending on s (respectively 
on t). Note that cos ip q — — - has two solutions (p q ^ and tp q ^ with | sin<^ g l | = | sm(p q ^\ = 



2tt 



d<p q 



[6(-sw<p qtl ■ (cp q - <p q> i))+5(-sin(p qi2 ■ (<p q - ^, 2 ))] 





2tt 



dp q S((p q -(p q>1 ) + 5(ip q -ip q>2 ) 



x* 



Thus 



sds dip q 



a 



VT 



a 



\t + ax\ 



(2.10) 



x- 



sign(t + ax) dt 



s>a, cos (fi q 



to 

arctan 



(t + ax) 2 + (1 + a 2 )(l - x 2 ) 

ti 



t + ax 



(2.11) 



to 



l + <? 2 ' v/(l + a 2 )(l-x 2 ) 
\t\ < (3 contains the upper limit of the t- integration as t < (3 =: t\. What concerns its lower 
limit t Q , the relation t > —f3 competes with t > —ax, as it follows from Eq. ()2.7|) for the 
lower limit a of the s-integration. This means for the ranges of the t- and ^-integrations, one 
has in the (vertical) stripe q = • • • l,p = 1 • • • oo of the g-p-plane to distinguish between 
the two regions, cf. Fig. 3: 

(a) region A with qp > 1 or a < (3 or —ax > —a > —(3, hence to — —otx, and 

(b) region B with qp < 1 or a > (3 or —a < —j3. 

In the case (b) one has again to distinguish between 

(i) x = — (3 /a ••• + /?/«, with the consequence to = —ax and 

(ii) x = (3 /a ■ ■ ■ 1, with the consequence t = —(3. 
Thus it is 



+ 1 +P 

dx J dt 

- i to 



+1 +p 
J dx J dt 

— 1 —ax 



for qp > 1 or a < (3, 



(2.12) 



+/3/a +0 +1 +P 

j dx J dt + J dx J dt for qp < 1 or a > (3. 

-P/a -ax +P/a ~P 



With the abbreviations (|2.5|) and 



f{t,x) 



a 2 — x 2 



arctan • 



t + ax 



v/(l + a 2 )(l-x 2 ) 
Eq. ()2.6j) can be written as (note that f(—ax,x) = 0) 



With 



1 oo 



X 1 = — dq dp 



16 7T 



1 



(p 2 -q 2 ) ( 



+ 6(1 - qp) 



+1 



+/9/Q +1 

dxf(j3,x)+ [ dx\f(j3,x)-f(-p,x)] 



-p/a 



-P/a 



+P/a 



+1 



+P/a 



dx{-l)f(-/3,x)= J dxf(P,-x)= I dxf(P,x) 

+P/a 
+ 1 



the terms for qp < 1 can be comprised as J dx f((3,x). Therefore 



(2.13) 



(2.14) 



(2.15) 



1 oo +1 

X 1 = - I dq I dp I dx -^-^ f(P,x) 



(2.16) 



1-1 

The final substitution p = 1/k transforms the region of the last two integrations from the 
(vertical) stripe g = 0- --l,p=l---ooto the more simple unit square q — • • • 1, k — • • • 1. 
With the abbrevations 

■ 2 1 - k 2 . „ 1 + q 2 k 2 



1 - q z 
2q 



2k 



> 0, a 



2qk 



> 1 



it is: 



i i +i 

f f 16 7T 1 

dq dk dx — f(p,x) w -30.70598 



(2.17) 



(2.18) 



oo-i 

The coefficients a, 8, and a of Eq. ()2.17|) make the integrand of ()2.18|) functions of q and k. 



Mathematica5.2 



17j yields the given figure. It seems to hold 



X X = -7T 4 



iln2-5# 

3 7T 2 



-30.705985239248893 • • • 



(2.19) 



How to derive this analytically ? Is this possible with the method of ref. 



III. THE INTEGRAL X 2 

The whole procedure of Sec. II is repeated here step by step. In Eq. (jA.ll|) . new variables 
<?' = (<?i + <?2)/ 2 > P' = (Qi ~ ?2)/2, and s' = e + (q 1 + q 2 )/2 lead to 

q x = q> + p', (e + Ql f = (s' + p'f, (e + q 2 f = (s' - p'f, 

q 2 = q'- p', (e + q 1 + q 2 f = (s' + q'f, e 2 = (s' - q'f. (3.1) 

Therefore 

2 3 1 

3„03 / z - 1 



X 2 = / cfg'dV 



g' 2 -p' 2 {q' + p') 2 {q' - p') 2 



x 



l ^ 5{{s> - q'f - l)e[(s> + q'f - 1]9[1 - + " - P') 2 ] • (3-2) 



The next step scales q', p', s' with A = 1 j\J\ — s' 2 according to q = Xq', p = \p\ s = As' 
with the consequences 

1 d 2 H 
l-s' 2 = ^-^, d 2 s'~ 



1 + s 2 ' (1 + s 2 ) 2 ' 

5((s'-q'f-l) = (l + s 2 ) 5(q 2 -2sq-l), 

±2sp + p 2 <l, 2sq + q 2 >l, -2sq + q 2 = 1, (3.3) 

from which follow q > 1 and p < 1 (what makes the energy denominator q 2 — p 2 positive) 
and sq > 0, s > a . Thus 

oo 1 +1 

16 7r f dx 8 q 2 p 2 



X2 = I dq I dp ^p-2 J Y (p 2 + q 2 ) 2 - (2qp) 2 X 



10 -1 

72, 



S — 5(q 2 - 1 - 2sq)Q(l - p 2 - 2sp)0(l - p 2 + 2sp). (3.4) 

1 + s 2 

s>a, cosip q >0 

Here and in the following of Sec. Ill the abbreviations 

q 2 — 1 - 1 — p 2 q 2 + p 2 

a = ^~: >0, (3 = —*->0, a = - - F > 1, 

2g 2p 2qp 

- = cos (f p = cos <(s,p), x = cosy? = cos<(q,p), y = cos </? g = cos <X(qr, s) (3.5) 
are used: 

oo 1 +1 

10 —1 s>o, cosip 9 >0 



t is introduced to replace s after having done the (p q integration. To this purpose the 
relation between t and s is needed. It follows from <p q + (p p = (p or (p p = (p — ip q , therefore 
cos(p p = cos(y? — ip q ) = cosy? cos <p q + sin 9? sin <£> g , what is in terms of t/s, x, and y (which 
equals +a/s because of the delta function): 



± y/l^sjl- or (t-ax) 2 = (I - x 2 )(s 2 - a 2 ) . (3.7) 



t a i - I (a 

- = x— 

s s 



This has the consequences 



S ~ a = 1_ X 2 » ( 3 ' 8 ) 

sds (t — ax)dt , . 

(3.9) 



1 + s 2 (t - ax) 2 + (1 + ct 2 )(l - x 2 ) ' 

With the help of Eq. ()3.8|) the <p q integration yields a function depending on s (respectively 
on t). Note that cos ip q = +j has two solutions ip q> i and ip qt2 with | sin<^ g l | = | sm<p q ^\ = 



1 - 



'a\2 



2n 



J -y [(K+sinv^i • {<p q ~ <p q ,i)) + 6(+smip qi2 ■ (ip q - (p q , 2 ))] 



o 

2tt 



dcp q 5((p q - cp q>1 ) + 5((p q - (p qj2 ) 2 v / T 3 ' 



;j; z 



o 



'o^2 



\/s 2 - a 2 



(3.10) 



Thus 



i 

sds dip q 1 /a \ yl — ^ 2 /" sign(t — ax) (it 



1 /tt \ v 1 — X 2 /" 

1 + s 2 2^ 6 V7 ~ V 0(/? ~~ = q J (t - ax) 2 + (1 + a 2 )(l - x 2 ) 

s>a, cosi/3 ? >0 to 



2 t — ax 

arctan ■ 



l + q 2 v/(l + a 2 )(l-x 2 ) 



(3.11) 

to 



\t\ < (3 contains the upper limit of the t- integration as t < (3 =: t\. What concerns its lower 
limit to, the relation t > —f3 competes with t > ax, as it follows from Eq. ()3.7|) for the lower 
limit a of the s-integration. This means for the ranges of the t- and x-integrations, one has 
in the (horizontal) stripe q = 1 • • • oo,p = • • • 1 of the g-p-plane to distinguish between the 
two regions, cf. Fig. 4: 

(a) region A with qp > 1 or a > (3 or —a < —(3, and 

(b) region B with qp < 1 or a < (3 or ax < a < (3, hence t Q = ax. 



In the case (a) one has again to distinguish between 

(i) x — — 1 • P/a, with the consequence t = —/3 and 

(ii) x = —/3/a ■ ■ ■ + /3/a, with the consequence t = ax. 
Thus it is 



+1 +p 
dx dt = < 



-1 to 



-/3/a +p +0/a +/3 

J dx J dt + f dx J dt for qp > 1 or a > [3 

-1 -f3 -0/a ax 

+1 +0 

J dx J dt for qp < 1 or a < (3 . 

— 1 ax 



With the abbreviations (jH.5j) and with 

2 



f(t,x) 



a 2 — x 2 



arctan ■ 



t — ax 



V /(1+«2)(1_ X 2 ) 

Eq. ()3.6|) can be written as (note that f(ax,x) = 0) 



X 2 

+ 9(gp- 1) 
With the identity 



00 1 

dq J dp 

1 o 

-P/a 



16 7T 1 



+ 1 



q 2 — p 2 1 + q 2 



e(l-gp) / da;/(^x) 



+y3/a 



dx[/(Ax)-M^)]+ / 



-/3/a 



-jS/a 



+ 1 



dx(-l)/(- / 3,x)= / dxm-x) 



dx f(/3,x) 



-i 



+/3/S 



+ 1 



the terms for qp > 1 can be comprised as J dx f((3,x). Therefore 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



dq J dp J dx 

I o -1 



16 7T 



q 2 — p 2 1 + q 2 



(3.16) 



This is similar to Eq. ()2.16|) . but there are also differences. The substitution q = 1/k 
transforms the region of the last two integrations from the (horizontal) stripe q — 1 ■ • • oo, p — 
• • ■ 1 to the more simple unit square k = ■ ■ - l,p = ■ ■ ■ 1. With the abbrevations 



a 



1-k 2 - 1-p 2 _ l + k' 2 p 2 
p = — z , a 



2k 



2p 



2kp 



(3.17) 



it is 



1 1 +1 

f f 16 7T k 2 

dk I dp j dx - - ^ ^—^ f{fl,x) . 



0-1 

Changing finally the notation with k —>■ q and p — > k makes a = a, (3 = j3, and a = a, cf. 
Eq. (|2.17jl . With these identities and with f(t, —x) = f(t, x) a further rewriting yields 
i i +i 



Xo 



dq dk dx 



16 7T 



(f 



1 - q 2 k 2 l+q' 



f(P,x) « 21.284906 



(3.18) 



oo-i 

This integral differs from Eq. ()2.18l) 'only' in an additional factor of — q 2 in the nominator 



of the integrand. Mathematica5.2 

"2 



Xo = 7i q 



■In 2 - 2 



171 ] yields the given figure. It seems to hold 
C(3)" 



7T Z 



21.284905670516334- •• . 



(3.19) 



How to derive this analytically ? Is this possible with the method of ref 



'• i 18 ] ? 



IV. THE CALCULATION OF X 

With Eqs. (pTSjl . (fTTT^ . (j3~TRjl . and with a, a, /3 being defined in Eq. (j2~T7jl the result for 
X = X x + X 2 is 

/",/"„ y , 16 7T 1 - g 2 2 /? + ax 

X = - dq dk dx — ±77 — arctan . = ss -9.42108 ■ • • . 

' ' ' 1 - q 2 k 2 l + q 2 a 2 -x 2 y/(l + a 2 )(l-x 2 ) 



0-1 



It seems to hold LJ, 

X = -7T 4 

The final result 



2, C(3) 
3 n z 



(4.1) 



-9.421079568732553 ■ • • . (4.2) 



S 2x = -^X = e 2x with -^- = iln2-^« 0.0241792 (4.3) 

appears in the weak-correlation limit of the Hugenholtz-van Hove (Luttinger-Ward) theorem 
|lJ,Ql- Because of /Z2x = e 2x it holds the sum rule fi 2x = S 2x analogous to fi x = E x . Whether 
perhaps also the more general expression T, 2x (k,uj) can be calculated in a similar way, has 
to be studied. 
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APPENDIX A: DERIVATION OF S 2x 

The one-body Green's function of the non-interacting system (ideal Fermi gas) 

(with O(x) = Heaviside step function) and G(k,u), the one-body Green's function of the 
fully interacting system, define the self-energy E(k, uS) through 

G(k, u) = G (k, oj) + G {k, w)E(Jfe, u)G(k, u) . (A.2) 



uj) appears in the Hugenholtz- van r 



ove theorem (in the Luttinger- Ward form \i — \iq 



„) with „ = chemieal potential, Q and in the GalitsMi-Migdal f Mm nla Q 

v = - [ d(k) 3 ! —^ 5 G{k,u)T,{k^) , 5 > . (A.3) 

2 J J 27T1 

v is the potential component of e, the total energy per particle. The contour of the 
(^-integration is to be closed in the upper complex tu-plane. In lowest order it is 
S x (fc) = —(1 + ^f~l n lTr|D- ^his ma kes v x = — f 2 ^ 1 , in agreement with v x = e x , what 



follows from the virial theorem v = r^-f-e. 



From the Feynman diagram for the exchange term of the 2nd-order self-energy it follows 

'"' r d 3 q x d 3 q 2 f dr\ x dr\ 2 



S 2x (fc, to) 



[atr. 



x 



4vr 4 J qfq'i J (2m) 2 
x G (\k + q 2 \,LO + r] 2 )G (\k + q x + q 2 \,to + r] 1 + r] 2 )G {\k + q x \, u + % 



(A.4) 



Use of (TO) yields 



E2x{k,uj) 



or. 



4tt 4 



d 3 q 1 d 3 q 2 



e(\k + q x + q 2 \ - 1)0(1 -\k + ^1)6(1 - \k + q 2 \) 
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\k 2 + q x ■ q 2 - iS 



6(1 -\k + Ql + q 2 \)e(\k + qi \- l)6(|fc + q 2 \ - I) 



to - \k 2 + q x - q 2 + i6 



(A.5) 



One may check this expression by using it in the Galitskii-Migdal formula (|A.3j) . Its lhs is 
known from the virial theorem as v 2x = 2e 2x with e 2x = energy in second order of exchange, 
calculated by Onsager et al. ^J. Its rhs gives with Eqs. (|A.1J) and (|A.5j) 

3(ar. q ) 2 f f d 3 kd 3 qi d 3 q 2 0(1 - fc)6(l -\k + q 1 + q 2 \)Q(\k + q x \ - l)Q(\k + q 2 \ - 1) 



rhs 



(2tt) c 



q x q 2 + \5 



+ 



/ 



d 3 kd 3 q l d 3 q 2 0(k - l)e(|fe + q x + q 2 \ - 1)6(1 - \k + q x \ - 1)6(1 - |fe + q 2 \) 



Qi ■ (~9 2 ) + iS 



(A.6) 



It is easy to show with the help of the substitutions q x 
that the second term equals the first one. Thus 



Qv Q2 -9' 2 , k -> — (fc + qf'J 



(27T) 



<?1 ■ 92 



(A.7) 



x 6(1 - k)Q(l -\k + q x + q 2 \)e(\k + q x \ - l)0(|fe + g 2 | - 1) = 2e 2x . 



P means the Cauchy principle value. This is in agreement with the above mentioned 
relation. That e 2x of Eq. (|A.7|) really agrees with the integral calculated by Onsager et 
al. || follows from the substitutions fe — > k x , q x — > q, q 2 — > — (k x + k 2 + q). From Eq. 
(|A.5|) also follows n 2x (k), the second-order- in-exchange contribution to the momentum 
distribution. It is again easy to derive the well-known asymptotics n 2x (k — > oo) = — g-^ ^Sr • 



After this control of £ 2x (/c, oj), the formula for E 2x = Re £ 2x (l, 1/2) follows from Eq. ()A.5|) 
as £ 2x = — ^pr (Ai + X 2 ), where X lj2 mean the integrals 

X x = [ d\ d 3 q 2 — — -L 0[i _ ( e + 9l + q 2 ) 2 ]0[(e + qA 2 - l]0[(e + q 2 f - 1] , 
J Qi ' Q2 Qi Q2 

(A.8) 

X 2 = [ d\ d 3 q 2 — — -^0[(e + 9l + q 2 ) 2 - 1]6[1 - (e + qi ) 2 ]0[l - (c + <Z 2 ) 2 ] • 

(A.9) 

They contain q 1 ■ q 2 as the energy denominator. l/qf 2 arises from the Coulomb repulsion 
and the remainder is due to the Pauli principle, e is a unit vector. Note that q l ■ q 2 < for 
X±, because of (2e + q^) ■ q x > and (2e + q 2 ) ■ q 2 > in combination with (2e + q x ) • q 1 + 
(2e + qr 2 ) • q 2 + 2g x • q 2 < 0. This latter inequality enforces q 1 ■ q 2 < 0. Thus X x < 0. It 
follows similarly X 2 > 0. The integrals X 12 do not depend on e. Therefore application of 
O = J d 3 e/ (4-7t) 2 S(e 2 — 1) does not change them (notice 2 5(e 2 — 1) = 5(e — 1)). Following 
Onsager et al. e is resolved into its components perpendicular to the q l — q 2 -plane e±, 
and in the plane ey: O = J de±/2 J d 2 e\\/(2-K) 2 5(ejj + — 1). The integration over 
may be done immediately by means of a change in scale: q 1 = q 1 y/l — e 2 ± , q 2 = q 2 a/1 — e\, 
ey = en-v/l — e\. The results are (denoting ey as e and deleting also all the other tildes for 
simplicity) 



P 



1 f ^ 9 
?2 V 2vr 



Ai = / cf gi cf g 2 

x6[l - (c + flfx + ^) 2 ]©[(e + q,? ~ l]0[(c + <Z 2 ) 2 - 1] , (A.10) 

X 2 = / Ai A 2 — ^- -J-j / p- 2 5(e 2 - 1) x 

X0[(e + 9l + q 2 ) 2 - 1]0[1 - (e + qi ) 2 ]0[l - (e + q 2 f] . (A.ll) 



Whereas the two terms of S 2x (fc, u), cf. (1A.6J) each contributes e 2x to t> 2x (thus t> 2x = 2e 2x ) 
as shown above, cf. (|A.7j) . X x and X 2 contribute different values to S 2x (l, 1/2) as shown in 
Sees. II and III. 



FIG. 1: The Feynman diagram of e2 x , analytically calculated by Onsager et al. 



a. 




FIG. 2: The Feynman diagram of E2 X (fc,o>), (semi) analytically calculated in this paper. 
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case a > [3 case a < (3 



FIG. 3: The dashed area is the region of integration described by Eq. 1)2.12(1 . 




case a > P case a < (3 



FIG. 4: The dashed area is the region of integration described by Eq. (|3.12j) . 



